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Abstract 

Although relative errors can readily be calculated, the absolute astrometric accuracy of 
the source positions in the Chandra Source Catalog (CSC), Version 1.0, is a priori 
unknown. However, the cross-match with stellar objects from the Sloan Digital Sky 
Survey (SDSS) offers the opportunity to compare the apparent separations of the cross- 
matched pairs with the formally calculated errors. The analysis of these data allowed us 
to derive a value of 0.16" for the residual absolute astrometric error in CSC positions. 
This error will be added to the published position errors in the CSC from now on, 
starting with CSC, Version 1.1. 
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1 Introduction 

The source positions in the Chandra Source Catalog (CSC; Evans et al., 2010) are 
characterized by error ellipses (circles for Version 1 of the CSC) which are based on the 
spatial distribution of the photons in the individual source detections. In the case of 
multiple detections of the same source in different observations the error ellipse of that 
source is derived from the error ellipses associated with the individual detections. These 
error ellipses provide a good measure of the statistical uncertainty of the location of the 
source in the frame of the observation, but leave out a series of potential sources of 
error that are external to the observation: 

• The error in the mean aspect solution for the observation; clearly, the effect of 
this error will be diminished when multiple detections of the same source are 
combined. 

• The calibration of the geometry of the spacecraft, in particular the optical axes of 
the aspect camera and the High Resolution Mirror Assembly (HRMA). 

• The astrometric errors in the Guide Star Catalog; this should be very small. 

• The calibration of the geometry of the focal plane, its projection on the detectors, 
and the distortions therein. 

For all practical purposes, we shall combine these errors and call it an astrometric 
systematic error, even though not all of its components are truly systematic. The intent 
of this study is to derive the value of this compound quantity in order to add it to the 
CSC statistical position error so as to obtain a reliable absolute error for each of the 
CSC sources. 
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Using the Sloan Digital Sky Survey (SDSS; York et al., 2000) object catalog (Abazajian 
et al., 2009) we have the opportunity to compare the formal statistical errors with the 
measured separations of CSC-SDSS cross-match pairs. A statistical analysis of these 
data allows us to determine an accurate value for the astrometric systematic error. 

This study is part of a larger project characterizing the contents of the CSC that will be 
described in more detail in a forthcoming publication (Primini, et al., 2010). 

2. Cross-identification 

We use the probabilistic algorithm of Budavari & Szalay (2008) to cross-match the CSC 
with the Seventh Data Release of the Sloan Digital Sky Survey (SDSS DR7; Abazajian 
et al., 2009). Using Bayesian hypothesis testing, one can objectively determine the 
quality of an association, which depends only on the measured positions and the 
astrometric uncertainties of the given sources. The Bayes factor is computed for every 
possible candidate association using a constant o s = 0.1" uncertainty for the SDSS 
sources and a varying o x positional error for each Chandra detection that is determined 
from the 95% accuracy limit s . In the high-precision regime, the dimensionless Bayes 
factor is calculated as 



where the uncertainties, as well as the angular separation \pl s between the two sources, 
are measured in radians. Thresholding on the above formula is not equivalent to cutting 
on the angular separation because of the varying uncertainties in Chandra. This has 
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proven to be superior for cross-matching GALEX and SDSS sources (Heinis et al., 
2009) that exhibit a similar behavior. 

Next we assign probabilities to the candidates based on a uniform prior that is 
determined from the ensemble statistics of the matched catalog in a self-consistent way, 
as described in Budavari & Szalay (2008) and Heinis et al. (2009). One may think of the 
prior as the average probability that a given source pair represents a physical match 
and it is expressed as 

N. 

P = 



N X .N S 

where N x and N s are the number of sources from either catalog in the intersection of the 
coverage of the two catalogs and JV* the number of true cross-match pairs, all three 
scaled to the entire sky. 

The posterior for each association is then given by 



P — 

r xs 



1 + 



i -pr 1 



B XS P 



that is reported in the cross-match catalog for each associations in addition to the 
angular separations and the Bayes factors. We note that, for a uniform P prior, a 
threshold on the posterior translates directly into a Bayes factor cut. However, the 
interpretation of the probability is much more straightforward. Consistency requires that 
N* is equivalent to the sum of the posteriors P xs over all source pairs. Most catalogs in 
the subsequent analysis use a P xs > 0.9. The CXC-SDSS cross-match catalog (version 
1 .0) is an exception, probably because a number of probabilities are underestimated 
due to the missing astrometric error that is the subject of this paper. However, we will 
apply the requirement P xs > 0.9 for the purpose of this study. 
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3 Procedure 

The CSC-SDSS cross-match catalog contains 7989 objects that are classified as stars 
in the SDSS catalog. Since these sources are, by their nature, point-like we assume 
their optical and X-ray positions to be well-determined and coincident. We have further 
narrowed the sample down by requiring the match probability to be greater than 90%. 
The resultant sample contains 6310 CSC-SDSS object pairs which are uniquely 
associated with 9476 source detections in individual observations; these 9476 objects 
were used for this analysis. By using the combined (CSC-SDSS) spatial error estimate 
of each object pair as the independent variable and analyzing the statistical distribution 
of the measured separations, it is possible to derive the value of the missing absolute 
astrometric error in the CSC. The assumption here is that the astrometric error is 
relatively small compared to the CSC uncertainties, especially off-axis, and will 
therefore mainly affect the pairs with small combined errors. What makes it possible to 
separate the astrometric error from the statistical error is the fact that the former is a 
constant, while the CSC statistical error varies over a wide range, primarily as a function 
of off-axis angle. 

The separation is a single-axis radial measure and, in order to perform the analysis 
correctly, the positional uncertainties also need to be converted to a single-axis radial 
quantity. CSC provides the major and minor axes of an error ellipse, while the SDSS 
gives independent errors in RA and Dec, which are also assumed to represent an error 
ellipse. However, in version 1 of the CSC the error ellipses are not fully implemented, 
yet, and instead approximated by circles (i.e., equal major and minor axes). The error 
ellipse in the SDSS is also close to a circle and thus the fact that SDSS did not report 
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covariance for RA and Dec errors is not very important. This justifies the use of a single- 
axis radial combined error. We derive this radial error by adding the geometric means of 
the major and minor axes for CSC and SDSS in quadrature; in other words: the square 
root of the sum (CSC plus SDSS) of the products of major and minor axis. We want to 
be dealing with 1-o values and since the CSC error ellipses refer to a 95% confidence 
level, the CSC values are to be multiplied by 0.408539. 

To state this in a more exact fashion, we define the following quantities: 

£o : semi-major axis of CSC 95% confidence ellipse 

£-1 : semi-minor axis of CSC 95% confidence ellipse 

Ora : 1 -o error in RA for SDSS positions 

o Dec : 1 -O" error in Dec for SDSS positions 

o c : 1-o combined statistical radial position error for CSC-SDSS cross- 
matches 

o a : 1-o astrometric error 

<j' c : 1-o combined corrected statistical radial position, including astrometric 
error 

p : (radial) separation of CSC and SDSS positions for a cross-match pair: 
measured error 

p N (o) normalized sample error 
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X 2 : reduced x 2 
Then four of these quantities can be expressed as: 

o c = 70.1669041 .£(,.£! + a RA .a Dec 
a' c = V0.1669041 .e .e 1 + a RA .a Dec + o% 

Pn O) = - 

~ 2 _ %i Pn 2 
X ~ n-1 

In the following, o c (or a' c ) is the independent variable, p or p N the dependent variable. 
All values are in units of arcsecond. 



4 Analysis 

After sorting the data in increasing order of on we calculated j 2 (pn (o c )) for bins of, 
successively, 100, 200, 300, 400, 500, 500, 500, and 476 sources, and plotted the 
results against the mean value of o c for each bin. The result is shown in Fig. 1a. The 
values at o c > 0.25 are quite reasonable, but the steep rise below this point is indicative 
of an error component that is of the same order. We interpret this as caused by the 
missing astrometric error discussed in the Introduction. Our assertion is that, if the left 
hand part of the curve can be flattened out by adding a suitable value for o a in o' c and 
using that value in the calculation of Pn(oc) and x 2 (Pn (oc)). one has determined the 
astrometric error. A value of o a = 0.16" (±0.01) provides a good result as shown in Fig. 
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1c. For comparison, the same plot for values of o a = 0.15" and o a = 0.17" is presented in 
Figs. 1b and 1d, respectively. 

To verify the reliability of the result, we plotted the distribution of p N in three ranges of 
the independent variable: o c < 0.15", 0.25" < o c < 1.0", and 1.0" < o c (Fig. 2). We expect 
these to show Rayleigh distributions; they do, but the one in Fig. 2a is significantly 
shifted toward higher values, as is to be expected. When we make the same plots 
again, using o' c and p N (a' c ) instead (see Fig. 3), the distributions all match. The 
distribution of Pn {p' c ) for the entire sample is shown in Fig. 3d. 

In Fig. 4 we present, for the bin sizes from Fig. 1, the average estimated error p against 
the average off-axis angle (in minutes of arc), including the 0.16 arcsec systematic 
error. As expected, small errors are predominantly found at small off-axis angles, large 
ones at large angles. 

Finally, in Fig. 5 we present, for the bins from Fig. 1 , the average source separation p 
against the average estimated error o' c , including the 0.16 arcsec systematic error; the 
dashed black line represents the identity relation. The figure shows that p tracks o' c quite 
well. But the divergence at higher values of o' c indicates that the statistical errors of the 
CSC positions are likely to be overestimated when those errors are large; this 
corresponds (cf. Fig. 4) to off-axis angles greater than 7 or 8 arcmin. The same 
phenomenon can be observed in Fig. 1c where the plot slopes down for large values of 
the error. 

We performed one more check on the results by calculating the x 2 function for varying 
ranges of off-axis angle 0. It appears that the function exhibits differences in slope on 
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the low-error side, depending on 0. In Fig. 6 we show the equivalent of Fig. 1 for all 
source pairs (dashed line) and for source pairs with 9 < 5D (solid line) assuming an 
absolute astrometric error o a = 0.18". However, drawing definite conclusions on the 
basis of this analysis is uncertain as yet, as the true nature of potential variations in 
positional errors across the field of view, and the associated change in PSF, is not fully 
understood. Consequently, there is not sufficient reason to change our recommendation 
of adopting o a = 0.16". 

5 Conclusion 

Our conclusion is that the astrometric error in CSC positions, resulting from the four 
components listed in the Introduction, is 0.16" ±0.02". Adding this value in quadrature to 
the statistical error associated with each individual detection will result in a reliable value 
for the absolute position errors in the CSC. This will be effected for the published 
position errors in all releases of the Chandra Source catalog, starting with Version 1.1. 
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Figure captions 

Fig. 1 x 2 (p N (cic)) and x 2 (Pn {o'c)) as a function of o c , respectively <j' c , for bins of 100, 
200, 300, 400, 500, ... , 500, 476 sources, (a) x 2 (p N (o c )) for uncorrected errors o c . (b) 
X 2 (p N (a c )) for corrected errors a' c , where o a =0.15 arcsec. (c) as panel b, for o a =0.16 
arcsec. (d) as panel b for o a =0.17 arcsec. 

Fig. 2 The distribution of p N (o c ) as a function of o c . (a) for o c < 0.15"; note that the 
Rayleigh distribution is clearly shifted to the right, (b) for 0.25" < o c < 1 .0". (c) for 1" < o c . 

Fig. 3 The distribution of p N {<j' c ) as a function of o' c for o a = 0.16". Note that all 
Rayleigh distributions agree within reasonable limits, (a) for o c < 0.15". (b) for 0.25" < o c 
< 1 .0". (c) for 1 " < o c . (d) for all values of o c 

Fig. 4 The measured position error (i.e., CSC-SDSS separation p; in arcsec) as a 
function of off-axis angle (in arcmin) averaged in bins of 100, 200, 300, 400, 500,..., 
500, 476 source pairs, ordered by 0. 

Fig. 5 The average value of the source separation, p, versus the average value of the 
estimated error (including a 0.16" astrometric systematic error) from the bins in Fig. 1 . 
The dashed black line indicates the identity function. The divergence at higher values in 
this plot (as well as the corresponding slope in Fig. 1c) hints that the position errors at 
off-axis angle greater than 7-8' (see Fig. 4) may be overestimated. 

Fig. 6 x 2 (pn {o' c )) as a function of o' c , as in Fig. 1, using o a = 0.18". The dashed line 
represents all sources and is similar to what is presented in the panels of Fig. 1 . The 
solid line represents only sources within an off-axis radius = 5D. 
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CSC-SDSS Crossmatch 
Reduced x 2 of Normalized Measured Position Errors 
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Normalized Error Distribution 
Uncorrected Errors 
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Normalized Error Distribution 
0.16" Systematic Error Added 
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